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ABSTRACT 

We present high resolution (R = 90,000) mid-infrared spectra of M dwarfs. The mid 
infrared region of the spectra of cool low mass stars contain pure rotational water 
vapour transitions that may provide us with a new methodology in the determina- 
tion of the effective temperatures for low mass stars. We identify and assign water 
transitions in these spectra and determine how sensitive each pure rotational water 
transition is to small (25 K) changes in effective temperature. We find that, of the 
36 confirmed and assigned pure rotational water transitions, at least 10 should be 
sensitive enough to be used as temperature indicators. 

Key words: Stars: low mass - stars: atmospheres - stars: fundamental parameters - 
infrared: stars 



1 INTRODUCTION 

Low mass stars (LMS) constitute ~80% of our stellar neighbourhood. They provide a probe of 
our understanding of main sequence stellar evolution and are key in determining the boundary be- 
tween stellar and sub stellar objects. However, their spectra are extremely rich in structure and their 
opacity is made up of many molecular and atomic absorbers, each with hundred of thousands to 
millions of spectral lines. This means that colours are not easily interpretable as diagnostic of their 
properties. Among the fundamental properties of LMS, of particular importance is their effective 
temperature, Teff'- the latter has frequently been investigated but it is still not well determined. In 
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particular there is not yet a tight correlation between spectral type and effective temperature as 
there has been a long standing discrepancy between empirical effective temperatures and those de- 
rived by synthetic spectra (~ 125 K as estimated by Kirkpatrick et al. 1993). This is primarily due 
to the lack of a complete inclusion of molecular and atomic opacities in the models, in particular 
in the near infrared (NIR). Modelling the atmosphere of LMS in the NIR is not a trivial task: water 
vapour dominates this part of the spectrum and at the effective temperature applicable to cool star 
atmospheres (<4000 K) water can access energies as high as 45000 cm~^ before it dissociates; 
to reproduce high temperature water spectra is very challenging because of the complexity of the 
vibrational and rotational motion of asymmetric triatomic molecules. Leggett et al. (2000, 2001) 
used a synthetic grid which includes the water linelist calculated by Partridge & Schwenke (1997) 
to determine fundamental parameters for a large sample of LMS. They conclude that problems 
remain with the match of the observed water bands. AUard et al. (2000) also find that the models 
including this water linelist reproduce well late-type dwarfs but fail to reproduce hotter dwarfs. 
The main sources of error in present water opacity databases are the incorrect high temperature 
transitions (Jones et al. 2002, 2005). 

Recently a new water linelist, BT2, has been produced (Barber et al. 2006). The BT2 linelist is 
by far the most complete linelist available, reaching higher energies than previously published. A 
comparison between a previous release of the BT2 linelist (BTl) and the PS linelist showed that at 
higher temperatures (T > 2000K) the PS linelist is missing around 25 per cent of the water vapour 
opacity (Jones et al. 2005). While the BT2 increases the accuracy of the fitting in the NIR region, 
its cut off at energies above 30000 cm~^ still means that none of the energy levels above 20000 
cm~^ are capable of being excited by 1 jira photons (Barber et al. 2006). 

Although the strongest water vapour opacity is at shorter wavelengths, analysis of sunspots 
(Wallace et al. 1995) shows that H2O also accounts for the majority of lines in the mid infrared 
(MIR) region between 8 and 21 jim. Polyansky et al (1997a,b) successfully assigned many of 
the water transitions in the 10-13 micron spectra of sunspots with further recent assignments be- 
ing made in conjunction with the analysis of oxy-acetylene emission spectra (Coheur et al. 2005; 
Zobov et al. 2006). In this spectral region the water transitions are largely pure rotational, in con- 
trast to the vibration-rotation transitions which dominate at shorter wavelengths. The use of pure 
rotational transitions greatly simplifies the spectral analysis as estimates of the transition strengths 
are much more straightforward. To a good approximation the strength of individual rotational tran- 
sitions can be estimated from a simple algebraic formula times a Boltzmann factor. Furthermore, 
the presence of pure rotational transitions within different vibrational states, which is characteristic 
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of the sunspot spectrum, yields a large dynamic range for temperature analysis. These advantages 
potentially make the MIR region ideal for the determination of the effective temperature scale for 
LMS. 

Here we propose a new methodology for establishing the effective temperatures of cool stars, 
in particular LMS, based on high resolution MIR spectroscopy. Taking advantage of the relatively 
easy to interpret MIR spectrum in cool stars has been done before: Ryde et al. (2006) used high 
resolution observations of water vapour in super-giants to show that classical photosphere models 
can not fit the observed spectra and that synthetic spectra based on cooler photospheric tempera- 
ture structure are needed. This supports the idea that individual rotational water lines are highly 
sensitive to temperature variations in this region. In this paper we present \hQ first high resolution 
(90,000) MIR observational data for two M dwarfs, together with one M giant. Our observations 
are described in Section 2. In Section 3 we present the data, the identification of the water lines 
and a theoretical sensitivity study using the BT2 linelist. In Section 4 we briefly conclude. 

2 OBSERVATIONS 

The data were taken in November 2006 using the Texas Echelon-cross-Echelle spectrograph (TEXES, 
Lacy et al. 2002) at the Gemini North 8m telescope. The observations were made on the nights of 
Nov 20, 23, and 27 for a total of ~13 hours. We observed three objects (BS587, GJ41 1 and GJ273 
- see Table 1) in the TEXES high spectral resolution mode in two grating positions: 1 1.2 and 12.5 
jira (800 and 892.857 cm~^). The two grating positions were chosen on the basis of a relative lack 
of narrow telluric features and abundance of temperature sensitive water features based on an early 
version of our ab initio water line list, as well as of a careful analysis of the sunspot (T~3300K) 
10-13 ^m spectrum. The spectral coverage was 0.5% at 11.2 fim and 0.75% at 12.5 jum. At 12.5 
yum, there are slight gaps in the spectral coverage because the angular width of the spectral or- 
ders is larger than the detector. The spectral resolution, as determined by Gaussian fits to telluric 
atmospheric lines in the 12.5 jxra setting, has a FWHM equivalent to R=90,000. We expect the 
same spectral resolution at 11.2 jira. The slit width was 0.5" for all observations and the length 
was roughly 4" and 2.5" at 11.2 jxra and 12.5 jira, respectively. For all observations, we nodded 
the source along the slit roughly every 10 seconds. 

Of our targets, GJ411 and BS587 were sufficiently bright to fine tune pointing based on the 
signal recorded from every nod pair. GJ273 was weaker and required accumulation of several 
nod pairs to evaluate pointing adjustments with confidence. On Nov 23, clouds were occasionally 
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Table 1. Properties of objects. 

Object RA Dec Sp. d (pc) v^„d (km/s) Mf,o! T<.//(K) logg 

GJ4111 11 03 20.2 +35 58 11.5 M2V 2.52 -84.74 8.88 3510 5 

BS5872 02 00 26.8 -08 31 25.9 M4.6III 182 +12.70 3425 -1 

GJ273-^ 07:27:24.5 +05:13:32.5 M3.5V 3.76 +18.10 9.56 3150 5 

^ Effective temperature taken from Tsuji et al. (1996); surface gravity from Jones et al. (1996); Mi,o; from Gushing et al. (2005). 
Effective temperature taken from van Belle et al., (1999); surface gravity from Fluks et al., (1994). 

Effective temperature taken from Tsuji et al. (1996); surface gravity from Jones et al. (1996); Mj,o( from Delfosse et al. (1998). 

present based on elevated background signal through TEXES. We discarded the data that were 
most severely effected. 

Coordinates and properties of the three objects are listed in Table 1. The two dwarfs were 
chosen among a sample of well studied dwarfs: note that for both GJ411 and GJ273, the quoted 
temperatures have an uncertainties of at least lOOK (cf Kirkpatrick et al 1993 and Tsuji et al. 
1996). The quoted radial velocities are from the literature (Marcy et al. 1987) and they match the 
sunspot features within 0.5 km/s. 

The data were reduced according to standard procedures as described in Lacy et al (2002). The 
TEXES pipeline provides for differencing of individual nod pairs, correcting spikes and cosmic 
rays, establishing a wavelength scale based on a user-identified telluric line, correcting distortions 
in the instrument, extracting ID spectra using optimal weighting along the slit length, and combin- 
ing spectra from the same night weighted by the square of their signal-to-noise ratio. The pipeline 
also performs flat-fielding and a first order correction for telluric spectral features using the differ- 
ence between an ambient temperature blackbody and sky emission. Further flat-field and telluric 
correction comes from division by a featureless point source such as an asteroid or hot star. Data 
from separate nights were aligned spectrally and combined according to the square of their signal- 
to-noise ratio using custom IDL procedures. The dispersion calculations are accurate to 0.6 km/s 
or better. 

3 ANALYSIS 

Figures 1 and 2 show the recorded spectra for our sample for the two wavelength regions. The 
gaps in the spectra coincide with the removals of bad pixels and of telluric lines. The signal to 
noise varies from star to star (~ 100 for BS587, 40 for GJ411 and 5 for GJ273). We used the 
bright star BS587 as a template for the identification of the water lines in the spectra due to its 
high S/N. In total we identified 71 lines in the observed spectra which could potentially be water 
(see Table 2); however, only 52 of them are confirmed by comparison with the BT2 linelist (see 
Table 2 and Figures 3 and 4), some of these (24, 25, 29, 38, 44, 53, 62 and 66 in Table 2) are weak 
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according to theory. Note that in order for a feature to be confirmed and assigned as water we 
required its intensity to be > 10^'^ cm/molecule and its frequency to be within 0.001 cm~^. 40 of 
the 52 Unes have been independently identified and assigned in the Sunspot spectra (Zobov et al. 
2006). Of the 19 remaining lines not identified in BT2, 13 are also present in the sunspot spectra 
(Wallace et al. 1995) and 4 of these 13 have been in fact assigned to water (Zobov et al. 2006 and 
private communications). Many of the observed features are blends of multiple lines (often 2, the 
para:ortho pair lines): for these cases we assigned the band to the strongest component (ortho in 
the case of pairs of lines). 

Figures 3 and 4 show examples of our line assignment with the observed stellar spectrum on 
the top panel overplotted to the sunspot spectrum (offset from the BS587 spectrum by +0.05), and 
with a synthetic water linelist computed at 3500 K on the bottom panel. Note that the choice of 
temperature for the synthetic water spectrum is not a best fit but was chosen to represent a bright 
cool stellar atmosphere and is close to that of the sunspot spectrum. In fact, a synthetic spectrum 
at 3000K may be more representative of a 3500K giant since the gas temperatures of the layers 
where most flux is emitted are typically 500K less than the effective temperature of the star (AUard 
et al. 1997). We have compared a 3000K synthetic water linelist with BS587 and find that the same 
number of water lines are identified. 

Some of the lines identified as water transitions in the observed spectra are not from experimentally- 
known levels. For some of these cases the quantum numbers have not been fully determined and 
those given in Table 2 are the rigorous ones only from the ab initio data. We have singled out these 
transitions in italic in Table 2. Finally, it is worth noting that some of the identified lines are in fact 
ro-vibrational: we have excluded these transitions from our sensitivity analysis (see next section 
and Table 3). 



3.1 Temperature sensitive water lines: a theoretical study 

Once the water lines have been assigned, we attempt to determine the sensitivity of each pure 
rotational transition to small (25 K) changes in temperature. This can be determined by simply 
estimating the relative intensity of the lines as a function of temperature by the use of modified 
Boltzmann equation: 

T -, ri— 15 Slower / Eiq^^^ 
Irel ~ 10 ""—^expi ^) 

where giower and Eiower are respectively the statistical weight and the energy level for the lower 
rotational state, and U is the total partition function (which is a function of the temperature, T). 
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Table 2. Transitions identified in the observed spectra of BS587. The confirmed water transitions are also given with energy levels, and quantum 
assignments when known experimentally. Transitions in italic are experimentally unknown transitions. The 10th column (Sunspot) states whether 
the confirmed water transition has been independentiy identified in simspots (Y) or not (N). The last column is the temperature of the maximum 
intensity gradient within the 1500 to 4000 K limits of our chosen model grid, (see Section 4). 
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Table 2. -continued 



Frequency (cm l) v\V2V3{\x) J (u) Ka Kc (u) 1^11^21^3 (i) J (1) Ka Kc (1) E(l) Sunspot Tmax 
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Note that our choice of 25 K is rather arbitrary: in order to use pure rotational water lines as good 
tools for the determination of the effective temperature of low mass stars, we require these lines to 
be sensitive to temperature changes below the temperature discrepancy found by other methods, 
usually at least lOOK. We calculated the partition function using the fitting formula of Harris et al. 
(1998) (Equation 12) which is suitable for high temperatures. We have considered temperatures 
ranging from 1500 K to 4000 K in steps of 25 K. Examples of intensities and gradients versus 
temperature are given in Figure 5 for two pure rotational water transitions: while the weak water 
transition at 797.424 cm~^ (line 5 from Table 2) exhibits a fairly low and flat gradient, the water 
transition at 894.637 cm~^ (line 58 from Table 2) is very sensitive to temperature changes up to 
~ 2800 K but it becomes less sensitive at higher temperatures. The temperature of the maximum 
gradient (i.e the temperature at which the transition is most sensitive to temperature changes) for 
all lines is also given in Table 2 (last column). We find that about 10 transitions out of the 36 
pure rotational ones identified in the BT2 linelist are sensitive enough to small (25K) temperature 
changes that they could potentially be used in the determination of an effective temperature scale 
for low mass stars. Clearly, a 25K sensitivity is a strict requirement. A 50K sensitivity would still 
make our methodology more accurate than others and would lead to the use of a larger number of 
lines: from the 36 pure rotational ones identified here we find about 14. 

The use of the most temperature sensitive water transitions for the determination of an accurate 
effective temperature scale of low mass stars becomes practical once a range of calibration objects 
has been observed with suitable measurements of water transitions. Our two dwarf stars only 
span around 300 K and have inadequate S/N for these purposes. Nevertheless, we attempt here 
to demonstrate the methodology by computing the intensity and gradient of the pure rotational 
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transitions from Table 2 at 3500K. Our computations are shown in Table 3 which includes the 
subset of rotational lines and identifiable in the Sunspot atlas. From this table we see that some 
of the most sensitive lines, at this particular temperature, are around 894.6 cm~^. In Figure 6 we 
plot the ratio of the two stellar spectra for a small wavelength region centred at 894.6 cm~^ and 
compare it to the ratio of two synthetic spectra computed at 3200 and 3500K. The strength of the 
'absorption' for some of the lines in the divided spectra is a clear indication that the transitions 
identified from Table 3 are indeed the most sensitive to temperature variations. It should be noted 
that the strongest most sensitive water lines occur towards lower temperatures. Observations of 
cooler objects should be within the reach of instruments such as EXES on SOFIA. The method 
can be calibrated with measurements of the growing number of cool low-mass eclipsing binaries. 

Finally, an effective temperature scale for low mass stars can only be quantitatively derived if 
effects such as stellar rotation, internal and atmospheric structure and equation of state are taken 
into consideration when generating the synthetic spectra. 



4 CONCLUSIONS 

In this paper we present the first high resolution MIR spectra of M dwarfs. Very high resolution 
observations of pure rotational water vapour transitions in the MIR may provide us with a new 
methodology in the determination of the effective temperatures for low mass stars. We have used 
the latest state-of-the-art water linelist (BT2, Barber et al. 2006) to identify and assign water tran- 
sitions in these spectra. In total we assign 52 water lines out of 71 likely water transitions; 36 are 
pure rotational lines. We have computed a theoretical sensitivity study to determine how sensitive 
each pure rotational water transition is to small (25 K) changes in effective temperature and we 
find that at least 10 should be sensitive enough to be used as temperature indicators. 

Due to the small wavelength region observed, the lack of a statistical sample of objects span- 
ning lower temperatures and the low S/N, it is impossible at this stage to determine the effective 
temperature of our two M dwarfs using the methodology outlined in this work; our technique is 
viable only provided we have a large enough number of sensitive rotational lines (hence a larger 
number of wavelength regions in the 8-21 ^m spectrum), as well as a large sample of objects span- 
ning the cool spectral subsequences. Both TEXES on Gemini and future generation telescopes 
such as the ELT have the potential to achieve this. 

Beside the need of a larger sample of objects, ultimately, for a realistic determination of an 
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Table 3. Relative intensity of pure rotational and assigned water lines at 3500K (Column 3), and Gradient at 3500K (Column 4). 
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effective temperature scale, model atmospheres including the BT2 water linelist should be used 
for future analysis. 
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Figure 1. Observed spectra centred at 11.2 fim. Top: BS587; Middle: GJ273; Bottom: GJ411. Spectra have been corrected for stellar velocities (see 
text). 




Figure 2. As for figure 1 for the 12.5 /nm region 
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Figure 3. Selected portions of spectra in the 11.2 jim region. Top: BS587 (Continuous) overplotted to the Sunspot spectrum (dashed, flux offset by 
+0.05 for clarity); Bottom: synthetic BT2 linelist computed at 3500K. 
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Figure 4. As in Figure 3 but for the 12.5 fj,m region 



14 S. Viti et al. 



797.424 cm-1 




—\ 1 r- 




_I I L_ 



1600 8000 8400 8800 3800 3600 4000 
Temperature in K 



1600 8000 8400 2600 3800 3600 4000 
Temperature in K 



894.637 cm-1 




1.2 




1600 8000 8400 8800 3800 3600 4000 
Temperature in K 



1600 8000 



8400 8800 3800 
Temperature in K 



3600 4000 



Figure 5. Relative Intensity (e.g. Column 2 in Table 2) versus temperature (Left) and actual gradient versus temperature (Right) for a weak (Top) 
and a strong (Bottom) pure rotational water transitions from Table 3. The top gradient has been multiplied by a factor of 10"^ to fit within the same 
scale as the bottom one. 
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Figure 6. Ratio of the GJ273 and GJ41 1 spectra (Top) compared to a ratio of two syntiietic spectra computed using the BT2 water hnehst, at 3200 
and 3500K. 
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